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Abstract 

We illustrate and emphasize the relevance of hyperbolic theories of 
dissipation in different physical scenarios. Particular attention is paid 
to self-gravitating systems where the relaxation time may become large 
enough as to require a description of the transient regime. It is argued that 
even outside that regime, hyperbolic theories may be needed to provide 
an accurate description of dissipative processes. 
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1 Introduction 



In the description of the evolution of any physical system, it is mandatory to 
evaluate, as accurateely as possible, the order of magnitude of different charac- 
teristic time scales, since their relationship with the time scale of observation 
(the time during which we assume our description of the system to be valid) 
will determine, along with the relevant equations, the evolution pattern. Take 
a forced damped harmonic oscillator and consider its motion on a time scale 
much larger than both the damping time and the period of the forced oscilla- 
tion. Then, what one observes is just a harmonic motion. Had we observed 
the system on a time scale of the order of (or smaller) than the damping time, 
the transient regime would have become apparent. This is rather general and 
of a very relevant interest when dealing with dissipative systems. It is our 
purpose here, by means of examples and arguments related to a wide class of 
phenomena, to emphasize the convenience of resorting to hyperbolic theories 
when dissipative processes, either outside the steady-state regime or when the 
observation time is of the order or shorther than some characteristic time of the 
system, are under consideration. Furthermore, as it will be mentioned below, 
transient phenomena may affect the way in which the system leaves the equi- 
librium, thereby affecting the future of the system even for time scales much 
larger than the relaxation time. 

2 The rationale behind hyperbolic theories 

Parabolic theories of dissipative phenomena have long and a venerable history 
and proved very useful especially in the steady-state regime [Q. They exhibit 
however some undesirable features, such as acausality (see e.g., Q, ||), that 
prompted the formulation of hyperbolic theories of dissipation to get rid of them 
0, Q. This was achieved at the price of extending the set of field variables 
by including the dissipative fluxes (heat current, non-equilibrium stresses and 
so on) at the same footing as the classical ones (energy densities, equilibrium 
pressures, etc), thereby giving rise to a set of more physically satisfactory (as 
they much better conform with experiments) but involved theories from the 
mathematical point of view. These theories have the additional advantage of 
being backed by statistical fluctuation theory, kinetic theory of gases (Grad's 
13-moment approximation), information theory and correlated random walks 
(at least in the version of Jou et al.) [[|. 

A key quantity in these theories is the relaxation time r of the correspond- 
ing dissipative process. This positive-definite quantity has a distinct physical 
meaning, namely the time taken by the system to return spontaneously to the 
steady state (whether of thermodynamic equilibrium or not) after it has been 
suddenly removed from it. It is, however, somehow connected to the mean col- 
lision time t c of the particles responsible for the dissipative process, ofentimes 
erroneously identified with it. In principle they are different since r is (con- 
ceptually and many times in practice) a macroscopic time, although in some 
instances it may correspond just to a few t c . No general formula linking r and 
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t c exists, their relationship depends in each case on the system under consid- 
eration. As mentioned above, it is therefore appropriate to interpret r as the 
time taken by the corresponding dissipative flow to relax to its steady value. 
Thus, it is well known that the classical Fourier law for heat current, 

q = - K VT, (1) 



with k the heat conductivity of the fluid, leads to a parabolic equation for tem- 
perature (diffusion equation), 

§ = XV 2 T ( X = JL) (2) 

at pc p 

(where x, P and c p are the diffusivity, density and specific heat at constant pres- 
sure, respectively), which does not forecast propagation of perturbations along 
characteristic causal light-cones (see j|, ||, [||, [0] and references therein). 
That is to say, perturbations propagate with infinite speed. This non-causal 
behavior is easily visualized by taking a look at the thermal conduction in an 
infinite one dimensional medium (see e.g. [0], [||]). Assuming that the temper- 
ature of the line is zero for t < 0, and putting a heat source at x = xq when 
t = 0, the temperature profile for t > is given by 



T 1 

1 oc — = exp 



(X - Xp 
t 



i2 



(3) 



implying that for t = =^ T = 5(x — xq), and for t = t > =^ T / Vx. 
In other words, the presence of a heat source at xq is instantaneously felt by 
all observers on the line, no matter how far away from xq they happen to 
be. The origin of this behavior can be traced to the parabolic character of 
Fourier's law, which implies that the heat flow starts (vanishes) simultaneously 
with the appearance (disappearance) of a temperature gradient. Although r 
is very small for phonon-electron, and phonon-phonon interaction at room 
temperature (0(lO~ n ) and 0(1O~ 13 ) seconds, respectively ||), neglecting it is 
the source of difficulties, and in some cases a bad approximation as for example 



in superfluid Helium |1C], and degenerate stars where thermal conduction is 



dominated by electrons -see Q, ||, [If |, for further examples. 

In order to overcome this problem Cattaneo and (independently) Vernotte 
by using the relaxation time approximation to Boltzmann equation for a simple 



gas derived a generalization of Fourier's law, namely 12] 



rf t +Q = ~^T. (4) 

This expression (known as Cattaneo- Vernotte's equation) leads to a hyperbolic 
equation for the temperature (telegraph equation) 

d 2 T dT 
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which describes the propagation of thermal signals with a finite speed 



v = Jx/r- (6) 



This diverges only if the unphysical asumption of setting r to zero is made. 

It is worth mentioning that a simple random walk analysis of transport 
processes naturally leads to telegraph equation, not to the diffusion equation 
-see e.g. [ 13 J . Again, the latter is obtained only if one neglects the second 
derivative term. 

It is instructive to write (0) in the equivalent integral form 





[ (t-t')l 


q = / exp 




T J-oo 


T 



VT(x, t')dt', 



(7) 



which in turn is a particular case of the more general expression 

q = - f Q{t-t')VT{x,t')dt'. (8) 

The physical meaning of the kernel Q(t — t') becomes obvious by observing that 

for Q = K5(t — t') =^ q = —kVT (Fourier) ^ 
for Q = constant =^ = x^ 2 T, 

i.e., Q describes the thermal memory of the material by assigning different 
weights to temperature gradients at different moments in the past. The Fourier 
law corresponds to a zero-memory material (the only relevant temperature 
gradient is the "last" one, i.e., the one simultaneous with the appearance of 
q). By contrast the infinite memory case (with Q = constant) leads to an 
undamped wave. Somewhere in the middle is the Cattaneo-Vernotte equation, 
for which all temperature gradients contribute to q, but their relevance goes 
down as we move to the past. 

From these comments it should be clear that different classes of dissipative 
systems may be described by different kernels. The one corresponding to (Q) 
being suitable for the description of a restricted subclass of phenomena. 

Obviously, when studying transient regimes, i.e., the evolution from a steady- 
state situation to a new one, r cannot be neglected. In fact, leaving aside 
that parabolic theories are necessarily non-causal, it is obvious that whenever 
the time scale of the problem under consideration becomes of the order of (or 
smaller) than the relaxation time, the latter cannot be ignored. It is common 
sense what is at stake here: neglecting the relaxation time ammounts -in this 
situation- to disregarding the whole problem under consideration. 
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3 Hyperbolic versus parabolic 



According to a basic assumption underlying the disposal of hyperbolic dissi- 
pative theories, dissipative processes with relaxation times comparable to the 
characteristic time of the system are out of the hydrodynamic regime. However, 
the concept of hydrodynamic regime involves the ratio between the mean free 
path of fluid particles and the characteristic length of the system. When this 
ratio is lower that unity, the fluid is within the hydrodynamic regime. When it 
is larger than unity, the regime becomes Knudsen's. In the latter case the fluid 
is no longer a continuum and even hyperbolic theories cease to be realiable. 

Therefore that assumption can be valid only if the particles making up 
the fluid are the same ones that transport the heat. However, this is (almost?) 
never the case. Specifically, for a neutron star, r is of the order of the scattering 
time between electrons (which carry the heat) but this fact is not an obstacle 
(no matter how large the mean free path of these electrons may be) to consider 
the neutron star as formed by a Fermi fluid of degenerate neutrons. The same 
is true for the second sound in supernuid Helium and solids, and for almost 
any ordinary fluid. In brief, the hydrodynamic regime refers to fluid particles 
that not necessarily (and as a matter of fact, almost never) transport the heat. 
Therefore large relaxation times (large mean free paths of particles involved 
in heat transport) does not imply a departure from the hydrodynamic regime 
(this fact has been streseed before [^4[, but it is usually overlooked). 

However, even in the case that the particles making up the fluid are respon- 
sible of the dissipative process, it is not "always" valid to take for granted that 
r and t c are of the same order (se e.g. [[Tq| , [|1|), or what comes to the same 
that the dimensionless quantity T = (tc s /L) 2 is negligible in all instances -here 
c s stands for the adiabatic speed of sound in the fluid under consideration and 
L the characteristic length of the system. That assumption would be right if r 
were always comparable to t c and L always "large" , but there are, however, im- 
portant situations in which r>t c , and L "small" although still large enough to 
justify a macroscopic description. For tiny semiconductor pieces of about 1CT 4 
cm in size, used in common electronic devices submitted to high electric fields, 
the above dimensionless combination (with r ~ 10~ 10 sec, c s ~ 10 7 cm/sec 
1 17]) can easily be of the order of unity. In ultrasound propagation as well as 
light-scattering experiments in gases and neutron-scattering in liquids the rele- 
vant length is no longer the system size, but the wavelenght A which is usually 
much smaller than L [18], [19]. Because of this, hyperbolic theories may bear 
some importance in the study of nanoparticles and quantum dots. Likewise 
in polymeric fluids relaxation times are related to the internal configurational 
degres of freedom and so much longer than t c (in fact they are in the range 
of the minutes), and c s ~ 10 5 cm/sec, thereby F ~ 0(1). In the degenerate 
core of aged stars the thermal relaxation time can be as high as 1 second [20|. 
Assuming the radius of the core of about 10 -2 times the solar radius, one has 
r ~ 0{Y) again. Fully ionized plasmas exhibit a collisionless regime (Vlasov 
regime) for which the parabolic hydrodynamics predicts a plasmon dispersion 
relation at variance with the microscopic results; the latter agree, however, with 
the hyperbolic hydrodynamic approach [21]. Think for instance of some syrup 
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fluid flowing under a imposed shear stress, and imagine that the shear is sud- 
denly switched off. This liquid will come to a halt only after a much longer 
time (r) than the collision time between its constituent particles has elapsed. 

The fact that r can quantitatively greatly differ from t c is most dramatically 
suggested by the matter-radiation decoupling in the early universe. In a recent 



paper Pavon and Sussman [22] by using the Lemaitre-Tolmann-Bondi metric 



1 23] along with the hyperbolic transport equation for the shear-stress 
1 + iTry 



rii cd h c X + n Q 



Tr] 



+ 2ij a ab = 0, with 77 = ^pr , (10) 



showed that the relaxation time of shear viscosity results several ordersof mag- 
nitude larger than the Thomson collision time between photons and electrons 
for most of the radiative era, i.e., in the temperature range 10 3 < T < 10 6 
Kelvin. (In above expression H a f, denotes the shear-stress tensor arising from 
the matter-radiation interaction, ij is the transport coefficient, and a ab the 
symmetric trace-free part of the gradient of the four-velocity, and p the ra- 
diation pressure. In this scenario t c is obtained by introducing the Thomson 
cross-section in Saha's formula. 

Many other examples could be added but we do not mean to be exhaustive. 

Even in the steady regime the descriptions offered by parabolic and hyper- 
bolic theories do not necessarily coincide. The differences between them in such 
a situation arise from (i) the presence of r in terms that couple the vorticity 
to the heat flux and shear stresses. These may be large even in steady states 
(e.g. rotating stars). There are also other acceleration coupling terms to bulk 
and shear stresses and heat flux. The coefficientsfor these vanish in parabolic 
theories, and they could be large even in the steady state, (ii) From the con- 
vective part of the time derivative (which are not negligible in the presence of 
large spatial gradients), (iii) From modifications in the equations of state due 
to the presence of dissipative fluxes ||. 

However, it is precisely before the establishment of the steady regime that 
both types of theories differ more importantly. It is well-known (see p], S, 



1 11], [|2j]) that a variety of physical processes take place on time scales of the 
order of (or even smaller) than the corresponding relaxation time, which as was 
stressed above does not imply that the system is out of hydrodynamic regime. 
Therefore if one wishes to study a dissipative process for times shorter than 
r, it is mandatory to resort to a hyperbolic theory which is a more accurate 
macroscopic approximation to the underlying kinetic description. 

In the particular case of neutron star matter, it appears that the relaxation 
time may indeed be of the order of magnitude of the characteristic time of some 
physical relevant processes. Thus, from (||) one has 

r = ~ 2 • 11 

v z pc p 

If the heat conductivity is dominated by electrons (as is the case of a neutron 
star), then we can adopt the expression [25] 
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K ~ 10 23 [ / o/10 14 g cm- 3 ][10 8 K/T]erg s" 1 cm" 1 K" 1 (12) 

and for the specific heat c p = (3T/M, where (3 (which is model dependent) is 
given by ]^6| /3 ~ 10 29 erg/K 2 and M denotes the total mass. Feeding back 



these expressions in (11) we get 

10 20 



[T 2 ] [v< 



(13) 



where [T] and [v] denote the numerical values of the temperature and the ve- 
locity of the thermal wave in Kelvin and cm/s, respectively. We have further 
assumed p ~ 10 14 g cm -3 and approximately 10 Km for the radius of the de- 
generate core. In this case, assuming for v the speed of light in vacuum, for 
temperatures of the order of [T] ~ 10 2 Kelvin -a certainly low value that corre- 
sponds to the latest phases in the evolution of a neutron star (see |26|)- we get 
r m 10~ 5 seconds. 

However, for more reasonable values of v, such as 10 3 cm/s, corresponding 
to the speed of second sound in superfluid helium, and for temperatures of 
the order of 10 9 Kelvin, one obtains r ~ 10 -4 seconds. Whereas for T ss 10 6 
Kelvin, the relaxation time may be found as large as r ~ 10 2 seconds. Therefore 
physical process on a time scale of the order (or smaller) than above belong to 
the transient regime and require a hyperbolic theory. 

Another astrophysical phenomenom where the effects of the transient regime 
can be felt is the double peaked temporal luminosity profile observed from x- 
ray bursters in star explosions, such as the one reported by Hoffman et al. p7| . 
These authors observed a x-ray burster showing a precursor peak, that lasted 
about 4 seconds, neatly separated from the main event -which lasted above 
1000 seconds. While the analysis of elastic photon diffusion through a cloud of 
ionized plasma of modest optical depth r* via parabolic theories just predicts 
the convential diffusion structure, the hyperbolic theory predicts a well defined 
peak for times shorter than tp ~ t 2 £ c /2 and a diffusive behavior for t > tjj 



This kind of event is reminiscent of laser induced heat pulses in solids where 



the double structure is evident [29]. Again hyperbolic theories successfully deal 
with such situations while parabolic do not Q. 

Finally, it is worth mentioning that relaxation time appears to be quite rele- 
vant in the outcome of gravitational collapse as many numerical and analytical 



calculations indicate (see [24], ]pO| l and references therein). In particular it has 
been shown that, for otherwise identical boundary and initial conditions, larger 
values of r imply more flattened and long lasting pulses of emission, thereby 
affecting the evolution of the object. 

4 Beyond the transient regime 

So far we have seen that only for times larger than r it is sensible to resort 
to a parabolic theory (modulo that the spatial gradients are not so large that 
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the convective part of the time derivative becomes important, and that the 
fluxes and coupling terms remain safely small). However, even in these cases, 
it should be kept in mind that the way a system leaves the equilibrium may be 
very sensitive to the relaxation time. 

Indeed, it has been shown [[Slf that after the fluid leaves the equilibrium, 
on a time scale of the order of relaxation time, the effective inertial mass den- 
sity of a dissipative fluid reduces by a factor that depends on the dissipative 
variables. By "effective inertial mass density" (EIMD) we mean the factor of 
proportionality between the applied three-force density and the corresponding 
proper acceleration (i.e., the three-acceleration measured in the instantaneous 
rest frame). The expression for the EIMD contains a contribution from dissipa- 
tive variables which reduces its value with respect to the equilibrium situation. 

Specifically, it has been shown that just after leaving the equilibrium on a 
time scale of the order of r the EIMD becomes (in relativistic units) 

A=(p + p)(l-a), (14) 



with 

kT 

a = ^Tp)> (15) 

giving rise, in principle, to the possibility of a vanishing A (i.e., a = 1) or even 
negative (a > 1). Such effect may be present in any dissipative fluid (either 
self-gravitating or not) j32fl . 

In order to evaluate a, let us turn back to coventional units. Assuming for 
simplicity p + p 2p, we obtain 

KT ^ M 42 

a = — — « — r—r- X 10 , (16) 

T(p+p) [t\[p\ 

where [k], [T], [t], [p] denote the numerical values of these quantities in c.g.s 
units. 

Obviously, except for extremely high values of k and T, a will be a much less 
than unity. Observe, in this connection, that although the smaller r, the larger 
a, this is of little consequence since ( |i~4|) is valid only on a time scale of the 
order of r. Therefore, a decreasing of A with physically relevant impact requires 
values of (|^) close to unity, due to large values of k and T but non-negligible 
values of r. 

A possible scenario where a may decrease substantially (for non-negligible 
values of t) might be provided by a pre-supernova event. At the last stages 
of massive star evolution, the decreasing of the opacity of the fluid, from very 
high values preventing the diffusion of photons and neutrinos (trapping [|33| ) , to 
smaller values, gives rise to a radiative heat conduction. Upon these conditions 
both k and T could become sufficiently large as to imply a substantial increase 
of a. Indeed, the values suggested in ||] ([k] ps 10 37 ; [T] ps 10 13 ; [r] « 10" 4 ; 
[p] « 10 12 , in c.g.s. units) lead to a ~ 0(1). The obvious consequence would 
be to enhance the efficiency of whatever expansion mechanism, of the central 



8 



core, because of the decreasing of A. Thus, it becomes clear that the value of 
a in the transient regime may critically affect the subsequent evolution of the 
system. Therefore it may be said that the future of the system at time scales 
much longer than the relaxation time (once the steady state is reached), may 
also critically depend on r. 



5 Concluding remarks 

It has been argued that the infinite speeds of propagation predicted by parabolic 
theories is just a consequence of using them outside their range of applicability 
which is limited by the discrete nature of matter In this connection it 

should be noted that hyperbolic theories also share this constraint but they 
naturally predict speeds in good agreement with experiments. Therefore one is 
led to conclude that hyperbolic theories have a range of aplicability much wider 
than parabolic theories. 

Aside from r hyperbolic theories introduce a certain number of new param- 
eters that couple the different dissipative fluxes. In spite of this being rather 
natural (as these theories are designed to explain more complex phenomena 
than parabolic theories do), this feature has been under unduly criticism as 
though these quantities were "free parameters" that one could choose at will 



1 15]. This is not the case. On the one hand, their number gets severely reduced 
when one realizes that they are interrelated. For instance, for an ideal gas un- 
der heat flux and bulk and shear stresses the parameters entering the transport 
equations of second order hyperbolic theory (equations (2.31)-(2.33) in [3. a]) 
are linked by the six equations (3.44)^(3.45) in [3. a]. On the other hand, these 
parameters are restricted by the convexity of the entropy function, and can be 
calculated explicitly with the help of kinetic theory, or fluctuation theory, such 
as in the case of a radiative fluid |3q] . 

In summary, although parabolic theories have proved very useful for many 
practical purposes, they appear to fail hopelessly in a number of well-known 
instances (such as transient regimes). By contrast, hyperbolic theories suc- 
cessfully predict the experimental results, and so, in these regimes hyperbolic 
theories happen to be more reliable. In the steady-state (under the conditions 
mentioned above) and for times exceeding r both theories converge. Before clos- 
ing we would like to emphasize we are not advocating to stop using parabolic 
transport equations. Our aim is simply to stress the convenience of using hy- 
perbolic transport equations there and when parabolic theories either fail or the 
problem under consideration happens to lie outside their range of applicability. 

We hope this paper will help to convince the reader that hyperbolic theories 
are indeed of not mere academic interest and it would not be wise to dispense 
of them. 
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